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ABSTRACT 


It  is  planned  to  send  an  infrared  scanner  in  a  balloon  up  to 
high  altitudes.  Azimuth  scanning  is  accomplished  by  the  torsion  of 
the  gondola.  The  elevation  angle  of  the  scanning  telescope  is  incre¬ 
mented  stepwise  after  each  full  revolution.  In  this  report  the 
scanning  mode  is  studied  and  improvements  thereof  are  suggested. 

For  a  follow-up  instrument  it  is  suggested  that  the  scanning  rate 
should  not  be  limited  by  the  slow  torsion.  In  a  next  approach  the 
e  evation  would  be  scanned  fast  compared  to  the  torsion.  Ultimately 
a  multiplex  scanner  should  look  at  many  directions  simultaneously. 
The  principle  of  a  multiplex  scanner  is  outlined. 
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Summary 


The  general  task  is  to  measure  from  a  balloon  the  infrared  radiation  S 
as  a  function  of  the  direction  it  comes  from,  as  a  function  of  time,  and  for 
our  (or  six)  bands  of  wave-numbers.  The  project  consists  of  three  parts: 

A.  Definition  of  the  optimum  scanning  mode  for  the  now  existing  LMSC  instrument 

B.  Design  considerations  for  a  future  instrument; 

C.  Development  of  the  software  for  the  evaluation  of  data  to  be  gathered  by 
the  existing  instrument. 

Part  A:  A  scanning  strategy  for  the  existing  instrument  is  outlined  in 
§§1-6.  Improvements  of  the  scanning  mode  are  suggested  in  §7.  Even  if 
feasible  most  of  these  suggestions  are  probably  too  late  for  thr  existing 
instruments.  Hence  the  corresponding  section  of  §7  belongs  actually  to  Part  B. 

Part  B:  In  the  long  range  the  rate  of  data  recording  should  not  be 
limited  by  the  slow  torsion  of  the  balloon,  which  provides  the  scan  movement 
for  the  azimuth  angle  a.  A  preliminary  solution  is  proposed  in  §7.1.4.  A 
long-range  attempt  for  increasing  the  data  recording  rate  is  based  on  the 
concept  of  "Fellgett's  multiplex  advantage".  The  principle  of  a  "multiplex 
scanner"  is  outlined  in  §9. 

Part  C:  The  aim  is  to  provide  a  mathematical  method  for  processing 
the  received  raw  data  to  convert  them  into  the  best  possible  "image"  of  the 
sky.  After  describing  the  data  flow  (§10,  11,  12)  a  linear  filter  approach 
for  image  restoration  is  presented  (§13).  Under  special  circumstances  the 
handicap  of  dilute  B-scanning"  (§12)  can  be  overcome  by  changing  the  receiver 
configuration  and  by  incorporating  a  nonlinear  scheme  into  the  digital  image 
restoration  procedure. 
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§1.  The  Objective 


The  goal  is  to  measure  infrared  radiation  (A  =  8-l2pm) .  xhe  measuring 
instruemtn  will  be  carried  by  a  balloon  up  to  a  height  of  about  30  km.  The 
instrument  consists  essentially  of  a  telescope,  a  point-receiver  in  a  ^xe 
position  of  the  image  plane,  and  a  recorder.  The  torsion  oscillation  of  t  e 
balloon  gondola  will  swing  the  telescope  around,  such  scanning  the  azimuth 
angle  a.  (For  details  about  nomenclature  see  §2.)  The  vertical  direction  of 
the  telescope  can  be  varied  under  remote  control,  such  pointing  at  different 
elevation  angles  8  on  the  sky.  A  filter  wheel,  also  under  remote  control, 
contains  various  color  and  neutral  density  filters. 


The  objective  is  to  devise  a  strategy  for  varying  the  elevation,  the 
filter  wheel  and  other  instrumental  parameters  which  can  be  optimized.  After 
establishing  the  nomenclature  (§2)  we  will  describe  the  relevant  features  of 
the  instrument  (§3).  Next  we  will  choose  a  suitable  coordinate  system  • 

Then  we  can  describe  the  scanning  strategy  (§5)  together  with  some  calibration 
and  housekeeping  procedures.  In  §6  we  will  discuss  the  choice  of  the  c  opper 
frequency.  In  §7  we  will  suggest  some  modifications  of  the  scanning  strategy, 
aimed  at  improving  the  efficiency  and  at  concentrating  on  the  most  interesting 
data.  §8  contains  some  systems  considerations,  and  §9  is  devoted  to  the 
"multiplex  scanner". 


The  rest  of  this  report  is  devoted  to  digital  image  restoration.  First 
we  describe  in  mathematical  terms  the  signal  flow  from  sky  to  tape  (slO).  e 
following  features  will  be  taken  into  account:  color  filters,  neutral  densi  y 
filters,  chopper,  lens  system,  receiver  size,  receiver  noise,  electrical 
filtering.  Next  (§11)  the  influence  of  the  neutral  density  filters  will  be 
considered  more  closely.  Since  the  step  AS  in  elevation  between  subsequent 
scanning  cycles  is  much  larger  than  the  resolution  68  the  recorded  signal  wil 
be  somewhat  ambiguous,  which  is  discussed  in  §12.  The  noise  which  is  added  o 
the  signal  in  the  receiver  can  be  suppressed  partially  by  subjecting  t  e 
recorded  data  to  a  digital  filtering  process,  which  is  a  generalization  of  the 
Wiener  filter  (§13).  Finally  we  consider  means  to  overcome  the  dilute  8 
scanning  (§12)  without  the  need  to  change  the  scanning  mode  of  the  first- 
generation  instrument.  This  method  will  work  if  all  object  details  (or  at 
least  most  of  the  interesting  ones)  are  points.  The  implementation  of  th 
scheme  requires  a  minor  modification  of  the  receiver  and  an  elaboration  of  the 
digital  image  restoration  method. 
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§2.  Nomenclature 


S(a,  3,  a,  t) :  radiation  coming  from  the  angular  intervals  (a,  a  +  da), 
(3,  3  +  dg)  in  the  spectral  interval  (a,  a  +  da)  during  the  time 
interval  (t ,  t  +  dt) . 

t:  time. 

A :  wavelength 

a  =  1/a :  wave  number. 

a:  azimuth  angle,  projected  into  horizontal  plane, 

a  ~  0:  magnetic  north  direction. 

1°  =  tt/180  rad  =  17.43  mrad;  1  mrad  =  0.0573  deg. 


ir/2  -  3  (in  rad):  elevation. 


3  =  0:  zenith. 


3  = 
a  = 
1/fi 

¥• 

w  - 


XM: 

6V 

L: 

v: 


V 


tt/2  rad:  equatorial  or  horizontal  plane. 

a„  ,sin{27rfi(t  -  t„)}:  torsion  of  gondola 
M  U 

:  full  period  of  torsion  oscillation. 

angular  amplitude  of  torsion;  maybe  >  2ir. 

2tt  a  ft  cos{2-rrfi(t  -  t  )}:  angular  velocity  of  gondola, 
M  0 

time  frequency, 
chopper  modulation  frequency, 
position  on  magnetic  tape, 
resolution  of  the  tape, 
total  length  of  the  tape, 
tape  velocity, 
bandwidth  of  magnetic  tape. 
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v^:  bandwidth  of  telemetry  system, 

index  E:  earth  related  coordinate  system. 


index  S:  sky  related  coordinate  system. 

index  G:  gondola  related  coordinate  system. 

3  (t) :  swiveling  of  the  gondola's  vertical  axis, 
s 

x',  y’:  orthogonal  coordinates  in  the  image  plane  of  the  telescope; 
the  x'-axis  is  horizontal. 

f:  focal  length  of  the  first  lens  in  the  telescope.. 

F  (cr) :  spectral  filter  functions,  m  =  1  .  .  ,.  4,  or  1  ,.  .  .  6. 

nx 

S  :  spectral  component  of  :S,  as  selected  by  E  . 

m  111 

z-axis:  gondola  axis  3  =  0* 
x-axis:  towards  magnetic  north, 

y-axis:  perpendicular  to  x-  and  z— axis . 

A3:  elevation  steps. 


63:  resolution  in  elevation 

T  :  transmittance  of  the  j-th  neutral  density  filter. 

3 

S  .(a,  3):  =  S  T  . 

mj  m  j 

M(x',  y’):  point  spread  function  of  the  lens  system,. 
M(v,  p) :  optical  transfer  function  of  the  lens  system, 
R(x',  y'):  spatial  distribution  of  the  receiver. 


P (a’si^^  3'):  angular  spread  function  of  the  telescope. 

S  (a) :  one-dimensional  equivalent  object  radiation. 

E 

G(v^) :  transfer  function  of  the  electrical  system. 

U  (t) :  signal  recorded  on  tape 
R 

D(vt):  digital  filter  function  for  noise  suppression. 


§3.  The  Parameters  of  the  Instrument 


?he  ?®lesc°Pe  ls  carried  Iby  a  balloon  to  an  altitude  of  about  30  km. 
where  it -.will  stay  for  about  four  hours.  The  elevation  of  the  telescope  nan  foe 
changed  in  flight.  In  vertical  position  (3  =  0)  the  telescope  would  point 
towards  the  balloon..  .  Hence  nothing  Is  lost  if  a  heat  source  for  a-callbratlon 
is  mounted  in  that  direction.  The  highest  elevation  (smallest  angle  -ft)  access- 
ible  for  radiation  measurements  of  the  sky  Is  called  30  •  The  Instrument  range 
f  3  extends  to  the  equatorial  plane  C3  =  tt/,2:;  also  called  parahorizon)  and  § 

Thl!!d'-f/he  ear.5h  h°riZon  at  about  %  =  r>2  +  1/10  (for  an  altitude  of  30  km), 
n  re  is  .danger  that  the  gondola  might  "swivel",  causing  3 -oscillation  of  the 

telescope  s  direction.  Nothing  definite  seems  to  be  known  about  the  amplitude 
and  frequency  of  these  3-oscillations.. 

Tihe  .azimuth  of  the  telescope  Is  f  ixed  with  respect  to  the  gondola.  Due 
f °  the  tors i°n  oscillation  of  the  gondola  the  telescope  will  aim  at 


a(t)  =  aM  sin{.2TTfif(t  -  ta)},. 


(3.1) 


The  torsion  amplitude  is  expected  to  foe  several  revolutions.,  say  aM (6tt.. 

It  can  be  predicted  only  roughly  before  the  flight..  Also  the  torsion  frequency 
J,£:S  irmerse>  tllle  torsion  period  1/13  cannot  foe  exactly  determined  from  pre- 
f light  measurements..  But  the  lower  bounds  for  aM  and  ft  should  be  known  before 
J  e  ignt.,  because  other  Instrumental  parameters  have  to  be  adjusted  accord- 

T°  ,8et  an  idea  °f  lt’he  orders  of  magnitude  we  start  from  the  assumption 
that  the  average  of  azimuthal  velocity  &  -  .o>  might  foe  ir  rad /min  -  1/20  rad/sec. 

hin.ce  ||  a)  |  AVERAGE  i*  we  get  with  aM  -  for  rad  now  1/ft  -  1500  sec  .24  min.. 

The  telescope  has  an  aperture  (diameter  of  5  Inches  .*»  125  mm,,  for  the 
longest  wavelength  <(12um)  this  results  in  a  theoretical  angular  resolution  of 
10  rad .  (Considering  the  difficulties  of  designing  the  infrared  lens  system 
consisting  of  so  few  elements  It  would  not  foe  surprising  to  measure  a  practical 
resolution  ten  times  larger,  say  6a  -  10" 3  rad  *  1  mrad..  Of  course  the  modu¬ 
lation  transfer  function  ((MTE)  would  .describe  more  completely  the  Influence  of 
the  lens  system.  In  the  Image  plane  of  the  telescope  are  four  ((not  only  one) 
receivers  horizontally  side-by-side.  Two  receivers  have  a  size  of  1x1  mrad2 
and  the  others  .2  x  .2  mrad2..  Their  angular  distances  are  each  5  mrad  from 
center  to  center,.  The  corresponding  length  measures  follow  from  the  focal 
ength  f  -  9, .5  Inch  240  mm  of  the  primary  lens .  In  the  Intermediate  Image 
plane  is  a  filter  wheel.,  which  contains  three  neutral  density  filters  of 
transmittances  10^,  1(0-2,  1(0-3^  and  four  color  fllter:S  spectral  ranges 

9-10,.5my ,  10,.5-11..5mp ,  and  11..5-12mp  .  Two  additional  color  filters 
(.8  ..3-9my ,  12-13my.)  are  available  as  standbys..  The  spectral  transmittance  of 
the  lens  system  and  the  spectral  sensitivity  of  the  receivers  are  known  to  be 
flat  within  1S.4~  The  temporal  response  of  the  receivers  is  characterized  by  a 
time  constant  of  2  msec..  The  noise  is  probably  caused  mainly  by  the  receiver 
independently  of  the  Incoming  radiation..  The  noise  .spectrum  decreases  somehow 
with  increasing  temporal  frequency  v,.  A  chopper  In  the  Intermediate  image 
plane  has  N  spikes  and  rotates  M  times  per  .second..  Hence  the  radiation  is 
modulated  by  yc  =  NM,  which  Is  planned  to  foe  .200  HZ.. 


The  signals  as  picked  up  by  the  four  receivers  are  recorded  on  a  multi¬ 
track  magnetic  tape,  which  moves  with  velocity  v.  The  resolution  length  fix 
and  the  total  length  L  of  the  tape  are  significant  parameters.  The  bandwidth 
of  the  tape  recorder  is  connected  with  fix  and  v  by  vM  =  v/6x.  The  dynamic 
range  of  the  tape  is  also  a  relevant  parameter.  Possibly  the  receiver  signals 
will  not  be  recorded  on  tape  but  transmitted  immediately  to  an  earth-based 
recorder.  In  that  case  the  bandwidth  v-p  of  the  telemetry  system  is  important. 


§4.  Choice  of  a  Coordinate  System 

I 

The  choice  depends  on  what  is  meaningful  for  describing  the  "object", 
what  dan  be  calibrated,  and  what  is  convenient  for  data  processing.  Basically 
one  could  choose  the  sky,  the  earth,  the  gondola  or  the  balloon  as  the  frame  of 
reference.  Since  the  earth  rotates  with  an  angular  velocity  of  2ir  rad/day  = 

4  mrad/min  the  stars  will  not  be  completely  stationary  objects  in  an  earth- 
related  system,  which  speaks  against  the  earth  frame  of  reference.  Nevertheless 
we  prefer  the  earth-related  coordinate  system  because  it  can  be  calibrated  most 
conveniently.  Gravity  provides  a  basis  for  definite  the  zenith  at  3  =  0,  and 
the  earth's  magnetic  field  can  serve  to  calibrate  the  azimuth  angle  a  =  0.  The 
earth  system  differs  from  the  gondola  system  due  to  the  torsion  and  due  to 
swiveling : 

aE  “  aG  =  aM  sin{2frfi(t  -  to)};  (4.1) 

“  ^G  ~  $s(t)  •  (4.2) 

The  center  (x '  -  0,  y'  =0)  of  the  telescopes  image  plane  corresponds  to  the 
point  a,  3  on  the  sky:  aj4  sin{2trfi(t  -  to));  3s(t).  A  shift  by  dx',  dy'  in  the 
image  plane  means  a  shift  in  angular  coordinates  by 

da  =  sin  g  dx'/f;  dg  =  dy'/f.  (4.3) 


The  various  spectral  filter  functions  are  called  Fm(a)  ;  m  =  1,  2,  3,  4. 
They  are  meant  to  describe  the  filter  transmittance,  but  also  the  lens  trans¬ 
mittance  and  the  receiver  sensitivity.  We  will  call  Sm  the  spectral  components 
of  the  object. 


Sm ( a  >  S,  t ) 


3,  a,  t)  Fm(a)  da. 


(4.4) 


§5.  Calibrations,  Scanning  Mode,  Housekeeping 


Before  the  flight  the  following  quantities  should  be  measured:  Trans¬ 
mittance  of  neutral  density  filters,  spectral  transmittance  of  neutral  density 
filters,  spectral  transmittance  of  color  filters  and  lens  system,  MTF  of  the 
lens  system  in  connection  with  various  color  filters,  flare  rejection  rate, 
spectral  sensitivity  of  receiver,  dark  current  of  receiver,  time-constant  of 
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receiver,  non-linear  response  of  receiver,  linearity  and  frequency  filter  of 
amplifiers,  linearity  and  frequency  response  of  tape,  chopper  frequency  vc . 

Ihe  parameters  a^,  ft,  gs,  which  relate  to  the  torsion  and  swiveling  of 
the  gondola  should  be  estimated  based  on  an  experiment.  A  method  for  choosing 
aM  deliberately  should  be  studied.  That  is,  the  gondola  might  be  pre-twisted 
and  hooked  before  take-off,  and  released  for  torsion  oscillation  when  at  opera¬ 
tional  altitude.  The  time  it  takes  to  perform  a  c-calibration  under  remote 
control  should  be  measured  in  the  laboratory. 

In  flight  the  instrument  should  be  checked  out  by  performing  at  first  a 
a-calibration.  For  this  purpose  the  telescope  has  to  be  pointed  at  the 
reference  heat  source  at  g  =  0.  Next  it  will  be  checked  if  the  azimuthal 
velocity  w  is  large  enough.  For  this  purpose  the  elevation  is  aimed  at  the 
earth  horizon,  where  many  radiation  sources  can  be  expected.  The  electrical 
signals  from  the  four  receivers  are  added.  The  four  signals  are  roughly  alike, 
but  delayed  by  5  mrad/o),  which  is  the  time  it  takes  for  an  object  point  to 
wander  from  one  receiver  to  the  next  one.  A  coarse  electronic  filter  bank  will 
find  out  the  resonance  frequency  to/ 5  mrad  in  the  combined  signal,  from  which 
a)  can  be  derived.  If  this  method  for  sensing  w  is  not  feasible  u)  should  be 
measured  by  some  other  method,  perhaps  based  on  centrifugal  forces.  At  any 
rate  it  would  be  a  waste  of  magnetic  tape,  if  the  measurement  would  begin  at  a 
very  low  torsion-angular  velocity  co. 

If  w  is  large  enough,  say  at  least  30  mrad/sec,  the  measurement  can  begin. 
The  telescope  is  elevated  to  the  steepest  angle  $o»  the  four  receiver  amplifiers 
are  disconnected;  the  tape  is  started;  the  starting  time  is  recorded;  the 
starting  angle  with  respect  to  a  magnetometer  Is  recorded  for  housekeeping 
purposes;  the  highest  transmittance  neutral  density  filter  is  inserted  together 
with  the  color  filter  m  =  1,  which  covers  the  whole  spectral  range.  During 
the  first  a-cycle  the  output  signal  of  the  amplifiers  is  checked  to  see  if  it 
exceeds  a  critical  level  which  corresponds  to  the  upper  end  of  the  dynamic 
range  of  the  tape  recorder.  If  this  level  was  exceeded,  the  next  denser  neutral 
density  filter  is  inserted  for  the  next  a-cycle.  Elevation  $o  and  color  filter 
F}  remain  the  same.  If  again  the  critical  level  is  exceeded  the  next  neutral 
density  filter  is  inserted  for  the  next  a-cycle.  If  now  an  a— cycle  is  completed 
without  exceeding  the  critical  level  color  filter  F2  is  inserted  and  the 
n.d.  filter  is  exchanged  for  the  next  higher  transmittance,  because  F^  lets 
pass  much  more  energy  than  any  of  the  other  filters  F2,  F3 ,  F4 .  The  latter 
three  are  roughly  alike  in  that  respect.  Once  F4  has  been  run  through  at 
elevation  $0  the  elevation  will  be  changed  $q  +  1/10  —  g^ .  Now  again  a  similar 
cycle  begins.  In  order  to  save  time  each  filter  F^  will  be  run  in  connection 
with  the  same  n.d.  filter  as  was  proper  at  elevation  g0.  This  guess  might  not 
always  be  adequate.  The  criterion  for  choosing  the  next  denser  filter  would 
be  as  before.  On  the  other  hand,  if  the  amplifier  output  never  comes  up  to  the 
critical  level  the  a-cycle  will  be  rerun  with  the  next  more  transparent  n.d. 
Filter.  For  "closing"  the  program  a  run  with  the  densest  of  all  n.d.  filters 
must  be  declared  as  successful  even  if  the  amplifier  did  exceed  the  critical 
level . 


The  increases  in  g-steps  by  1/10  rad  will  go  on  till  g  has  reached  an 
elevation  where  the  telescope  aims  at  the  earth.  That  g-position  will  not  be 
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executed;  rather,  the  command  will  be:  return  to  and  start  all  over.  A 
series  of  cycles  over  all  0-positions  takes  much  longer  time  than  a  half¬ 
period  1/2  ft  of  the  torsion  oscillation.  The  measurement  cycle  will  be  inter¬ 
rupted  when  the  angular  velocity  w  falls  below  critical  level.  During  the 
off-time  (while  | to f  <  ice)  the  a-calibration  is  repeated.  The  duration  of  the 
"off-time"  depends  on  ft  and  toc.  Giving  these  parameters,  or  at  least 

conservative  estimates,  allows  specifying  the  time  available  for  the  a-cali¬ 
bration  and  other  instrumental  check-ups. 

A  housekeeping  track  on  the  tape  shall  record  clock-signals,  a-raarkings 
from  the  magnetometer  and  records  of  all  events  like  changes  of  filters, 
elevation  jumps,  and  calibration  results. 


§6.  Choice  of  the  Chopper  Frequency 


It  is  advantageous  to  chop  the  radiation  periodically,  because  a.c. 
amplifiers  are  easier  to  operate  than  d.c.  amplifiers.  Also  the  signal-to- 
noise  ratios  might  be  influenced  favorably,  as  we  will  discuss  after  establish¬ 
ing  two  boundaries  for  the  chopper  frequency  vc.  The  chopper  frequency  vc  is 
the  product  of  the  revolutions  per  unit  time  of  the  chopper  wheel  and  the 
number  of  angular  periods  on  the  rotating  chopper  mask. 

Vc  has  to  be  high  enough  so  that  the  signal  is  still  recognizable  after 
chopping.  The  finest  periodic  object  structure  which  can  be  resolved  has  a 
period  of  6a  rad.  When  scanned  by  the  telescope  which  spins  with  angular 
velocity  oj  rad/sec  this  periodic  object  structure  will  generate  in  the  receiver 
a  periodic  signal  with  a  time  frequency  u)/6a.  The  chopper  frequency  vc  which 
acts  like  a  modulation  carrier  frequency  has  to  be  twice  as  high  as  the  highest 
ohject  frequency,  hence 

vc  >  2o)/5a.  (6.1) 

In  the  worst  case  oj  will  be  u  =  ^MAX  =  ^  •  Assuming  6a  =  10“3  rad,  a^j  = 

6 tt  rad,  1/ft  =  1500  sec,  we  get 

vc  >  4 it  aM  ft/ 6a  =  160  HZ.  (6.2) 

On  the  other  hand  the  chopper  frequency  should  not  be  unnecessarily 
large.  Otherwise  the  bandwidth  of  the  tape  might  not  be  sufficient.  The  upper 
sideband  of  the  modulated  signal  has  a  frequency  of  vc  +  iu/6a  (c.f.  equ.  6.1). 
The  highest  possible  frequency  vty  to  be  recorded  on  the  tape  depends  on  the 
tape  speed  v  and  the  tape  resolution  6x  like  =  v/6x.  Hence  we  must  request 

vc  +  w/6a  <  =  v/6x.  (6.3) 

Since  we  want  the  tape  speed  to  be  as  low  as  possible  we  select  the  smallest 
possible  chopper  frequency  vc  (equ.  6.1  and  6.2), 
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3io/6a  <  v/6x. 


(6.4) 


mmm 


Making  the  same  assumptions  as  for  equ.  (6.2),  and  also  6x  =  25m  it  follows 
as  minimal  tape  velocity 

v  >  3cow.v  <5x/<5a  =  6mm/sec.  (6.5) 

MAX 

Over  a  period  of  4  hours  =  14,400  seconds  this  requires  a  tape  length  of 

L  =  v  •  14,400  =  86.4  m.  (6.7) 

Now  let  us  discuss  how  the  signal-to-noise  situation  influences  the  choice  of 
the  chopper  frequency  vc .  Although  the  specific  assumptions  are  rather  vague 
this  consideration  is  nevertheless  meaningful  in  principle.  The  additive  noise 
spectrum  might  be  proportional  to  1/v,  and  the  sensitivity  of  the  receiver 
might  vary  like  1 / [ 1  +  (v/vR)2].  Hence  the  signal-to-noise  ratio  at  the 
chopper  frequency  is 

vc/[l  +  (vc/vR)2].  (6.8) 

It:  has  a  maximum  if  vc  is  chosen  to  be  equal  to  vr,  at  which  frequency  the 
receiver  sensitivity  has  dropped  to  1/2  of  its  d.c.  maximum.  Probably  v^, 
which  is  roughly  the  inverse  of  the  receiver  time  constant,  will  be  at  least  in 
the  KHZ  region,  so  it  is  not  sensible  to  bring  up  the  chopper  frequency  to  the 
same  level.  Yet  equ.  (6.8)  confirms  at  least  qualitatively  that  it  is  better 
to  chop  than  to  work  with  d.c.  signals. 


§7 .  Improved  Scanning  Mode 

The  scanning  mode  described  in  §5  is  based  on  the  discussions  with  AFCRL 
and  LMSC  representatives  (August  1969) .  That  scanning  mode  should  be  compatible 
with  the  hardware  design.  However,  several  aspects  could  be  improved  with  a 
moderate  amount  of  hardware  changes.  The  improvements  aim  mainly  at  better 
light  efficiency  and  matching  the  bandwidths  of  the  various  components. 


§7.1  (a,  6)  Scanning  Mode 

Employing  the  scanning  mode  as  described  in  §5,  it  takes  a  fairly  long 
time  to  complete  full  (a,  8)  cycle.  It  is  mainly  the  angular  velocity  u>  of 
the  gondola  which  slows  down  the  data  rate.  The  electronics  and  the  tape 
recorder  could  digest  much  higher  data  rates.  Avoiding  this  mismatch  and  intro¬ 
ducing  some  conveniences  are  the  aims  of  this  section. 


§7.1.1  B-Interlace 


The  8-steps  were  designed  to  be  0.1  rad,  while  the  8-tesolution  is 
perhaps  100  times  finer.  This  leaves  fairly  wide  8-zones  unobserved.  To  elim¬ 
inate  this  handicap  at  least  partially,  the  first  (a,  8)  cycle  might  cover 
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m  w 


rings  3o»  3o  ^  0*1*  3q  +  0.2,  etc.  tne  second  (a,  3)  cycle  would  be  interlaced 
like  3o  +  0.01,  3o  +  0.11,  3o  +  0.21,  etc.,  and  the  third  cycle  3g  +  0.02, 

30  +  0.12,  etc.  This  numerical  example  is  adequate,  if  total  time  (4  hours) 
would  cover  ten  (a,  3)  cycles. 


§7.1.2  Concentrating  on  the  Horizon  Belt 


The  most  interesting  details  are  probably  to  be  found  around  the  para- 
honzon  3  =  tt/2  ±  1/10.  Hence  within  this  belt  the  3  increment  should  be  finer 
than  at  the  higher  elevations. 


§7.1.3  Continuous  3-Increase 


At  the  moment  when  3  is  increased  by  a  discrete  amount  (like  0.1)  some 
data  might  be  lost,  because  the  telescope  might  vibrate  for  a  short  while. 

Hence  a  constant  d3/dt  might  be  better.  If  an  increment  of  33  per  revolution 
is  wanted  dB/dt  should  obey 

(d3/dt)/(da/dt)  =  63/ 2tt ;  da/dt  =  u).  (7.1) 

This  implies  that  m  -  w(t)  must  be  known.  Knowledge  of  m  can  be  gained  in 
several  ways.  t(the  clock)  and  a(t)  are  recorded  anyway  for  housekeeping 
purposes,  w  can  be  derived  from  oi  and  t.  Or  t  and  w(t)  are  constantly  measured. 
Or,  since  the  time  t  is  not  interesting  explicitly  it  would  be  good  enough  to 
monitor  a  and  u(a).  At  any  rate  (d3/dt)/(da/dt)  would  be  <<  1,  because  a  high 
ds/dt  of  the  telescope  would  probably  cause  mechanical  problems. 


§7.1.4  Vertical  Telescope  and  Rotating  Mirror 


So  far  the  data  rate  depends  mainly  on  the  angular  velocity  to  of  the 
torsion  and  on  the  angular  resolution  6a.  Pro  unit  time  the  telescope  will 
point  at  to/ 6a  different  resolution  elements.  With  to  =  tt  rad/min  =  1/20  rad/sec 
and  6a  =  1/1000  rad  this  results  in  50  data  per  second.  Such  a  low  data  rate 
is  certainly  far  below  what  can  be  handled  conveniently  by  the  receivers,  the 
electronics  and  the  tape  recorder.  This  data  rate  is  probably  also  far  below 
a  level  where  photon  noise  would  be  severe.  Hence  we  suggest  a  different  (a,  3) 
scanning  mode  with  a  much  higher  data  rate,  which  can  be  adjusted  freely. 

The  telescope  is  fixed  in  the  gondola.  It  points  vertically  to  the 
zenith  (3=0).  A  prismatic  body  with  quadratic  or  polygon  cross  section  has 
mirror  surfaces.  This  multi-mirror  is  placed  on  top  of  the  telescope  where  it 
rotates  around  its  axis  in  horizontal  position.  A  horizontal  hollow  cylinder 
in  front  cf  the  rotating  mirror  selects  coarsely  the  3-range  of  interest.  The 
3-scan  rate  d3/dt  is  two  times  the  angular  velocity  of  the  spinning  mirror. 

The  3-scan  rate  can  now  be  much  larger  than  the  a-scan  rate  a).  The  a-steps 
depend  on  the  ratio  oi  w  to  d3/dt,  and  of  the  number  of  facets  on  the  spinning 
mirror.  This  number  has  to  be  matched  to  the  8-belt  of  interest,  which  is 
preselected  by  the  hollow  cylinder.  The  four  receiver  elements  may  be  lined 
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up  in  parallel  to  the  axis  of  the  spinning  mirror.  Hence  they  scan  at  different 
a-positions,  avoiding  redundancy  of  the  four  receiver  signals  as  in  the  origianl 
setup.  Originally,  the  reason  for  letting  the  four  receivers  scan  through  the 
(a,  3)  domain  at  exactly  the  same  paths  was  the  possibility  for  detecting  fast 
temporal  changes  of  the  object  S(a,  3,  a,  t).  This  reason  is  not  very  valid 
anymore,  because  with  the  vertical  telescope  a  full  (a,  3)  cycle  will  be 
finished  much  faster,  perhaps  a  hundred  times. 


§7.1.5  Monitoring  the  Swiveling  of  the  Gondola 


Instruments  for  measuring  the  swiveling  3s(t)  of  the  gondola  are  probably 
available  in  the  aircraft  industry.  But  it  could  be  chat  3s(t)  is  difficult 
to  measure  because  3S  is  very  small  (hopefully)  and  very  slow.  In  that  case  a 
simple  hydrodynamic  instrument  might  be  desirable.  A  tube  system  could  form 
a  triangular  scar  as  large  as  the  gondola  floor.  The  ..three  spikes  are  slightly 
bent  upwards.  The  tube  consists  of  electrically  insulating  material  and  is 
partially  filled  with  a  conductive  fluid.  Pairs  of  electrodes  are  at  the  ends 
of  the  three  spikes.  Whenever  a  spike  tilts  down  the  fluid  will  close  the 
electrical  connection  in  that  spike.  Monitoring  the  on-off  events  might  be 
sufficient  enough  for  deducing  the  swiveling  of  the  gondola. 


§7.2  Modifications  of  the  Chopper  and  Filter  Wheel 


The  reasons  for  modifying  the  chopper  and  filter  wheel  are  two-fold: 
some  objects  on  the  sky  might  reveal  a  clue  about  their  identity  if  their 
spectral  characteristics  are  known,  for  example  a  quantity  like  (9S/3a)/S 
To  get  this  quantity  we  might  want  to  measure  the  spectral  components  AVERAGE 
Si,  S2 ,  S3 ,  S4  in  rapid  sequence  before  (a,  3)  have  changed  markedly.  Further¬ 
more  the  efficiency  of  collecting  radiative  energy  and  data  can  be  increased 
by  making  the  scanning  mode  more  flexible,  instead  of  designing  it  rigidly  on 
the  basis  of  worst-case  considerations. 

§7.2.1  Instant  Adaption  of  Neutral  Density  Filters 

So  far  (  5)  we  had  planned  to  change  neutral  density  filters  only  after 
full  a-cycles.  That  simple  strategy  might  require  that  the  same  a-cycle  be 
repeated  several  times  when  n.d.  filters  are  required.  This  waste  of  time  can 
be  avoided  if  the  n.d.  filters  change  instantly  whenever  certain  clipping 
levels  are  reached.  The  hardware  to  do  this  can  be  designed  as  in  a  modern 
recording  photometer. 


§7.2.2  High-Speed  a-Scanning 


Instead  of  exchanging  color  filters  stepwise  after  full  a-cycles  one 
might  place  little  pieces  of  color  filter  in  sector-form  on  the  chopper  wheel. 
Remember,  the  chopper  wheel  is  imaged  1:1  onto  the  receiver  elements.  For 
s implying  the  argument  let's  say  the  filters  used  so  far  have  each  a  spectral 
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passband  of  Imp  width,  covering  the  ranges  (8-9mp) ,  9-lOmp),  (lO-llmp)  and 
(H-12mM) .  It  would  be  better  to  take  a  set  of  complementary  filters  instead: 
(9-12mu),  (8-9  +  10-12mu) ,  (8-10  +  ll-12my) ,  (8-llmy) .  This  does  not  improve 
the  color  discrimination,  but  it  increases  the  average-color  signal  by  a  factor 
three,  which  is  desirable  if  the  scanning  speed  is  energy-limited. 


§7.2.3  Angular  Color  Wedge  Filters 

Color  wedge  filters  (OCLI)  are  based  on  the  fact  that  (in  good  approxima¬ 
tion)  the  cr-transmittance  curve  of  a  multilayer  film  will  be  stretched  along 
the  0-axis  by  a  factor  M,  if  all  layer  thicknesses  are  reduced  by  a  factor 
1/M.  Angular  color  wedges  are  produced  by  evaporating  onto  a  substrata  which 
spins  with  variable  velocity.  A  thin  sector  slit  exposes  only  a  small  portion 
at  a  time  for  the  deposition  of  a  layer.  Hence  thickness  and  cr  will  vary  with 
the  angular  position  cp  on  the  angular  color  filter. 

The  angular  color  wedge  might  be  easier  to  put  onto  the  chopper  wheel 
than  sector  pieces  of  interference  film.  Also  a  continuous  a-variation  is 
probably  more  convenient  than  stepwise  o-changes  for  data  evaluation.  It 
should  be  possible  also  to  produce  an  angular  color  wedge  with  a  "complementary 
passband"  as  described  in  §7.2.2. 


§7.2.4  Angular  Color  Wedge  for  Fourier  Spectrometry 


If  later  the  number  of  resolvable  spectral  elements  should  be  larger  than 
four  (as  now)  the  multiplex  advantage  of  Fourier  spectrometry  might  be  exploited 
in  the  following  way:  First  a  fixed  pre-filter  selects  the  total  range  (8-12mp). 
Next  an  angular  wedge  filter  with  transmittance  like  1  +  cos(b  o  cp)  covers  the 
receivers  and  spins  (cp  -+  cpt)  *  The  o-integrating  properties  of  the  receiver 
produce  the  Fourier  transform  of  the  source  S(a)  at  that  particular  point 
(a,  6) •  To  say  the  same  thing  in  another  way,  the  spectral  transmittance  of  the 
spinning  filter  is  such  that  each  spectral  component  a  is  chopped  by  its . own 
specific  chopper  frequency  vc(a).  Effectively  the  same  can  be  achieved  if  a 
linear  color  wedge  of  the  cosine-variety  is  placed  on  top  of  a  slit  receiver. 

A  sector  star  with  its  number  of  spikes  increasing  in  proportion  to  the  radium 
spins  on  top  of  the  slit  receiver  such  that  each  point  has  its  own  specific 

chopper  frequency. 


§7.3  Economy  Considerations  for  Receivers,  Amplifiers  and  Recorder 


As  pointed  out  earlier,  the  data  rate  in  the  receiver  is  so  low  that 
probably  the  bandwidths  of  the  four  receivers,  the  four  amplifiers  and  the  tape 
recorder  are  not  fully  utilized.  The  setup  with  the  vertical  telescope  (§7.1.4) 
was  already  a  step  in  the  right  direction  towards  more  economical  bandwidth 

matching . 
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§7.3.1  Only  One  Amplifier 


To  have  four  amj.  fliers  in  parallel  is.  not  only  four  times  as  expensive, 
it  also  requires  balancing  these  amplifiers,  for  otherwise  one  could  not  compare 
the  outputs  from  different  amplifiers.  Since  the  bandwidth  of  an  amplifier  can 
be  easily  much  wider  than  the  temporal  spectrum  produced  by  scanning  and 
chopping  of  the  signal  one  can  assign  different  carrier  frequencies  to  different 
object  parameters  (time  frequency  multiplexing).  One  example  of  this  category 
is  the  Fourier  spectrometer  (§7.2.4).  Other  examples  would  require  other 
chopper  disk  configurations. 


§7.3.2  Fresnel-Encoding  Receiver 


The  receiver  may  consist  of  a  narrow  long  stripe  which  is  covered  by  an 
absorbing  pattern  like  a  sector  of  a  Fresnel  zone  plate.  Since  there  is  only 
one  receiver  only  one  amplifier  is  needed.  The  scanning  direction  is  now  along 
the  receiver  axis  such  that  each  object  point  creates  a  temporal  "chirp"  or  FM 
signal.  This  recording  method  compares  to  a  point-receiver  method  like  Girard's 
grille  spectrometer  compares  to  a  slit  spectrometer.  The  decoding  of  the 
chirp-convolved  signal  can  be  done  by  a  digital  computer,  by  an  electronic 
correlator  or  by  an  optical  analog  computer,  which  is  nothing  but  a  holographic 
reconstruction  setup.  Remember,  a  Fresnel-hologram  is  the  convolution  of 
object  and  Fresnel  pattern.  Image  reconstruction  means  to  undo  the  Fresnel 
convolution.  That  is  what  we  have  to  do  with  a  photographically  recorded 
amplifier  output.  Since  this  artificial  hologram  has  been  convolved  with  a  one¬ 
dimensional  Fresnel  pattern,  a  lens  system  which  includes  at  least  one  cylinder 
lens  is  needed,  as  it  is  well  known  from  one-dimensional  holography. 


§7.3.3  Use  of  an  A/D  Converter 


An  A/D  converter  behind  the  amplifier  has  several  advantages.  In  this 
project  it  is  probably  most  important  that  non-linear  distortions  of  any 
analog  recording  medium  are  circumvented.  An  A/D  converter  would  also  be 
natural  if  the  data  shall  be  transmitted  in  PCM.  The  question  might  arise  if 
temporal  modulation  by  a  chopper  is  meaningful  if  an  A/D  converter  follows. 

The  answer  is  probably  "yes",  because  the  signal  frequencies  should  be  shifted 
away  from  the  noisy  frequency  region  around  d.c.  Finally  an  A/D  converter 
would  be  natural,  if  the  system  shall  communicate  with  a  digital  computer. 


§7.3.4  Variable  Speeds  of  Chopper  and  Tape  Recorder 


If  the  data  rate  is  proportional  to  angular  torsion  frequency  w,  the 
chopper  might  slow  down  if  the  torsion  gets  slower:  vc  *  w.  Since  chopping 
is  something  like  sampling,  it  means  that  the  chopper  sampling  rate  is  always 
sufficiently  high  to  match  the  data  scan  rate.  But  if  the  chpper  frequency 
varies  proportionally  to  w,  the  tape  velocity  might  do  the  same  in  order  not 
to  waste  any  tape.  A  tape  recorder  velocity  proportional  to  w  =  da/dt  would 
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be  convenient  because  now  the  length  coordinate  x  on  the  tape  is  directly 
proportional  to  the  angular  coordinate  a  of  the  object  S(.a;. 


§8  Systems  Considerations 


§8.1  Systems  Layout 

During  the  actual  balloon  experiment  various  types  of  data  will  flow 
between  the  various  components  as  depicted  in  figure  1.  The  a  in 

this  figure  are  explained  in  table  1,  "List  of  Components  .  There  are  the 
data  signals  gathered  by  the  LMSC  instrument  from  where  they  are  sent  to  an 
A/D  converter.  Then  there  are  several  command  signals  which  initiate  ac  ion 

like  removing  the  cap  from  the  telescope.  Furthermore  ^f^he  state  of  the 
signals,  which  record  angular  and  time  coordinates  as  well  as  the  stat 
instrument  (filter  wheel,  etc.).  The  purpose  of  devising  figure  l  was  to  get 
a  feel  for  the  total  system  and  to  make  sure  that  enough  housekeeping  signal 
are  recorded  so  that  the  data  signals  can  be  properly  identifi  . 


§8.2  Economy-Reliability  Considerations 


However,  a  suitable  recorder  must  be  on  board.  The  ne  of1|_he  onents 

are^othe^factor^which"  influence^  the  strategy  of  the 

experiment . 

There  are  nine  conceivable  options  depending  on  whether  a  control  unit 
is  up  and/or  down  and  whether  a  storage  unit  is  required  up  and/or^dow:  ^  ^ 

These  nine  options  together  with  the  equipment  requ  flexible. 

?  ThP  ninth  option  is  most  expensive  but  also  most  itiidu 

IfcRL  chose  rtis  option  in  order  to  be  best  prepared  for  many  unforseeable 
situations  during  the  first  experiment. 


§8.3  Test  Procedures 

Before  the  actual  experiment  th^ins trument  tas^to  t^tested^in^variou^ 

Si«iX  JjU-  and 

^a^^r^frr  t^cr :;dWUral  effects  such  as 
chopper  frequency  and  electronic  time  constants. 
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§8.4  Remaining  Problems  to  be  Covered  in  §§10-13. 

A  mathematical  formulation  of  the  flow  of  signal  data  has  to  be  developed 
taking  into  account  the  modulation  transfer  function  (MTF)  of  the  telescope  and 
the  frequency  responses  of  the  electronic  components.  Based  on  this  mathematical 
model  a  computer  program  has  to  be  devised  which  allows  recpvery  of  the  original 
data  by  removing  the  influence  of  the  MTF,  etc.  as  much  as  >is  meaningful  under 
the  signal-to-noise  circumstances. 


Table  1 

List  of  Components 


Cap 

CAP 

Elevation  drive  (3 1 

[ 

ELEV 

Chopper  | 

‘  LMSC  Telescope 

CHOP 

Filter  Wheel 

FILT 

Magnetome tome ter 

MAGN  a 

da/dt  (d/dt  of  magnetometer  or  centrifugal) 

da/dt 

Clock  t 

CLOCK 

Altimeter 

ALTI 

Telem.  Recv.  and  Araplif. 

TELE-R 

De-Multiplexer 

DE-MLPX 

On  Board  Control  Unit 

CONTROL 

Amflif.,  A/D 

A/D 

8-Track  Store 

8-TRACK 

Multiplexer 

MLPX 

Amplif.  and  Telem. 

Emitter 

TELE-E 

Swivel  Monitor  8S(+) 

SWIV 

4  Radiation  Receivers 

4  RAD  RCV 

D/A  or  A/D  Converter 

D/A  or  A/D 

Calibration  Source 

CAL 

Emitter 

EMIT 

D/A  and  Multiplexer 

D/A  &  MLPX 

Computer 

COMP 

Receiver 

RECEIV 

Store 

STORE 
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CAL  SWIV 


GROUND 


Table  2 

Nine  Instrumental  Options, 


FUNCTION 

Up 

Down 

Both 

EQUIPMENT 

Up 

Down 


§9  A  Multiplex  Scanner  Based  on  the  Talbot  Effect 


In  1836  F.  Talbot  described  an  effect  which  was  explained  about  fifty 
years  later  by  Lord  Rayleigh.  Today  this  eff  ect  is  known  by  the  names  ''Fourier 
images"  ( Cowley  and  Moodie,  Proc  Phys.  Soc.  B70,  486 ,  1957;  B7 6,  378,  1960) 
or  "self-imaging"  (Montgomery,  J.  Opt.  Soc.  Am.  57,  772,  1967).  This  effect 
has  been  used  for  Fourier  spectrometry  (Lohmann,  Proc.  I.C.O.  Conf,  on  Opt. 
Instr.,  London,  1961).  We  will  briefly  describe  this  Fourier  spectrometer  and 
then  explain  how  it  can  be  generalized  into  a  spectrometric  scanner.  By 
"spectrometric  scanner"  we  mean  an  instrument  which  collects  information  as  a 
1  jnct  ioii-'df  wavelength  X  and  of  the  source  coordinates. 

Consider  figure  °,  but  for  the  moment  not  quite  as  general  as  shown.  A 
monochromatic  poin*-  j.  rce  may  be  located  on-axis  (not  off-axis  as  in  figure  2) 
in  the  front  focal  plane  of  a  lens.  The  resulting  plane  wave  exp(i  X  Z)  falls 
onto  a  grating  in  plane  Z  =  0.  The  complex  light  amplitude  immediately  behind 
the  grating  is 

U(X..  0)  =  1  +  cos(27TX/d)  =  1+1/2  e+lKAX/d  +  i/2  e"lKXX/d.  (9.1) 

Herein  K  is  as  usual  K  =  2tt/A.  Th^  second  form  of  U(X,  0)  is  meaningful 
because  each  term  causes  a  plane  wave  in  the  half-space  Z  i  0  behind  the  grating 

U(X,  Z)  =  e1KZ  +  1/2  eiK(Xi/d  +  ZA-(X/d)2) 

+  1/2  e1K(-Xi/d  +  2/1-(i/d)2).  (9.2 ) 

The  validity  of  this  plane  wave  representation  of  the  wavefield  U(X,  Z)  can  be 


checked  easily  by  inserting  it  into  the  wave  equation  (V2  +  K2)U  =  0.  We 

re-write  U(X,  Z)  into  a  form  which  is,  convenient  for  computing  the  intensity 
I(X,  Z;  A)  =>  |U(X,  Z)  | 2 , 

U(x,  Z)  -  oIK7-  [1  +  cos  (2irX/d)  .-ilU<  1  -W?' I , . 

I(X,  Z;  A)  =  1  +  cos2(2?!X/d)  +  2  cos(2nX/d)  cos[KZfi  -  /T. .  j] 

=  3/2  +  1/2  cos  (4/rX/d)  +  2  cos(2uX/d)  cos[KZfl  -  /TTTjJ; 
cos[KZfl  -  /...)]  &  cos(irZA/d2) .  (9.3) 


The  important  feature  about  the  intensity  distribution  I (X ,  Z;  A)  is  its 
longitudinal  periodicity  with  a  period  AZ, 


AZ - ZL - 

k{  l  -  /7T7} 


(9.4) 


In  other  words  at  distances  AZ,  2AZ,  3AZ  the  intensity  is  the  same  as  in  Z  =  0 
immediately  behind  the  grating.  Hence  it  is  justified  to  speak  about  "self¬ 
imaging  of  a  grating  which  has  been  illuminated  by  a  monochromatic  plane  wave. 
At  certain  intermediate  planes  where  the  "modulation  factor"  cos (ttZA /d2)  is 
minus  one  the  contrast  of  i(X,  Z;  a)  will  be  reversed,  while  at  yet  other  planes 
the  modulation  factor  will  be  zero,  which  means  that  i(X,  Z;  A)  has  no  contrast 
in  the  fundamental  lateral  period  d. 

Now  let  us  consider  how  this  "Talbot  effect"  can  be  used  for  spectrometric 
purposes.  Suppose  the  light  source  consists  of  a  doublet,  A  and  A  +  i).  = 

A (1  +  1/N).  The  total  intensity  I(X,  Z;  A)  +  I(X,  Z;  A  +  6A )  will  consist  of 
two  parts  with  slightly  different  longitudinal  periods  AZ(A)  =  2d2/A  and 
AZ (A  +  <5A)  =  2d2 / ( A  +  6A).  Both  parts  coincide  close)y  behind  the  grating  at 
Z  =  0,  but  they  will  get  increasingly  out  of  step  at  larger  distances  Z.  In 
particular  at  Zq  =  Nd  /A  the  two  parts  will  cancel  each  others'  modulation  since 
the  two  modulation  factors  have  opposite  signs, 

cos  ( irZ^/d2)  =  cos(tiN); 

cos{7rZ0(A  +  6A)/d2}  =  cos  { 7T  (N  +  1)}.  (9.5) 


Hat ( e  the  doublet  separation  6A  =  A/N  can  be  deduced  from  the  distance  Zq  at 
which  zero-modulation  is  observed:  5A  =  d2/Z0. 


For  a  source  with  a  more  general  spectral  distribution  S ( A ) 
intensity  will  be 


I(X, 


Z;  A)  dA  =  I(X,  Z). 


the  total 


(9.6) 


Measuring  the  lateral  modulation  of  I(X,  Z)  for  various  distances  Z  yields 
essentially  „he  same  information  as  a  Fourier  spectrometer: 


cos(TiZA/d2)  dA. 


(9.7) 
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The  only  difference  is  that  the  Fourier  transform  is  carried  out  as  a  function 
of  A  and  not  1/A  as  usuel .  The  actual  measurement  of  the  modulation  of 
I(X,  Z)  is  done  by  placing  a  second  grating  at  a  distance  Z  behind  the  first 
grating  (see  figure  2).  The  second  grating,  which  has  the  same  grating  constant 
d  is  moved  laterally  with  velocity  vx  and  the  transmitted  intensity  is  collected 
and  measured  by  a  receiver  RCV.  The  amplitude  of  the  a.c.  part  of  the  photo¬ 
electric  signal  with  the  time-frequency  vx/d  represents  the  modulation  term 

j S (A)  cos ( ttZA /d 2 )  dA. 

Now  let  us  consider  how  such  an  instrument  can  gather  information  also 
about  the  geometrical  structure  of  the  source  S(a,  3,  A).  For  simplicity  let 
us  assume  for  the  moment  a  monochromatic  point  source  but  no  longer  on-axis. 

If  the  point  source  is  shifted  vertically  no  change  occurs  in  the  intensity 
distribution  I(X,  Z;  A)  behind  the  first  grating.  (This  statement  is  correct 
only  if  ZA3/d4  <<  1.)  But  if  the  point  source  is  shifted  laterally sby  an 
amount  f  tan  a  the  intensity  distribution  will  be  "skewed"  to  I(X  +  Z  tan  a, 

Z;  A).  Hence  the  horizontal  shift  would  show  up  as  a  phase  delay  in  the 
temporal  a.c.  signal. 

In  this  form  our  instrument  would  "sae"  only  the  horizontal  but  not  the 
vertical  structure  of  the  source.  In  order  to  overcome  this  limitation  we 
propose  to  spin  the  whole  instrument  around  the  Z-axis  with  an  angular  velocity 
to.  At  the  same  time  we  move  the  second  grating  laterally  with  velocity  vx  and 
also  longitudinally  with  velocity  vg  because  we  want  to  measure  the  x-modulation 
at  all  distances  Z.  It  is  plausible  that  three  types  of  movements  are  needed 
because  we  convert  the  three-dimensional  source  information  S(a,  6,  A)  into  a 
one-dimensional  time  signal,  the  photoelectrical  current.  Each  movement 
introduces  a  temporal  carrier  frequency  in  the  photoelectrical  signal.  We 
started  to  investigate  how  the  choice  of  the  movement  parameters  to,  vx,  vz  will 
ultimately  influence  the  resolution  in  three  dimensions  of  S(a,  6,  A).  The 
proper  choice  of  these  parameters  is  also  important  in  order  to  make  sure  that 
the  three-dimensional  object  S(a,  8,  A)  can  be  deduced  unambiguously  from  the 
one-dimensional  photoelectrical  current. 
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§10  Signal  Flow  from  Sky  to  Tape 


The  object  signal  is  the  radiation  S(a,  8,  a)  coming  from  the  sky.  S 
might  also  vary  as  a  function  of  time.  But  we  assume  these  temporal  variations 
to  be  slow  compared  to  the  azimuthal  scanning  cycle  of  the  instrument.  The 
m-th  color  filter  Fm(o)  selects  the  spectral  component  Sm 

Sm(a,  8)  =  f S(a,  8,  a)  Fm(a)  da.  (10.1) 

The  j-th  neutral  density  filter  reduces  the  signal  by  a  factor  T j ; 

Tj  Sm  (a  5  8 )  —  *-*mj  ^  >  8 )  •  ( 10  •  2) 

The  chopper  acts  multiplicatively  like 


Smj  (ot , 


27riHvrt 
e  c 


(10.3) 


Only  the  fundamental  chopper  frequency  ( | £ |  =  1)  will  later  pass  through  the 
electronic  amplifier.  Hence  it  is  enough  to  consider  the  | & |  =1  term  as  the 
signal  which  enters  the  optical  system  of  the  telescope, 

Smj(a,  8)  2  Aj  cos(2Trvct) .  (10.4) 

Now  we  must  relate  the  sky  coordinates  a,  8  to  the  coordinates  x',  y' 
in  the  image  plane  of  the  telescope.  An  auxilliary  coordinate  system  x,  y,  z 
is  defined  with  its  origin  in  the  center  of  the  telescope's  lens  (see  figures 
3,  4  and  5).  The  z-axis  is  also  the  gondola's  axis  pointing  towards  the 
zenith  8=0.  The  x-axis  points  to  the  magnetic  north  pole  at  azimuth  angle 
a  =  0,  and  the  y-axis  is  perpendicular  to  both  x-  and  z-axes.  Presently  the 
telescope  axis  may  be  pointed  towards  aj  ,  6i •  Hence  the  sky  coordinates 
ai ,  8i  are  projected  into  the  center  x'  =0,  y'  =  0  of  the  telescope's  image 
plane.  A  sky  point  close  to  aj ,  8i  in  a  =  aj  +  a',  8  =  8i  +  8'  will  occur  at 
x'  =  a'f  sinBj ,  y'  =  8'f  in  the  image  plane.  The  letter  f  stands  for  focal 
length.  We  have  used  the  small  angle  approximation  tan  a'saa',  tan  8 'f«8 *  »  which  is 
justified  since  the  receivers  in  the  image  plane  are  all  very  close  to  x'  =  0, 
y'  =  0.  The  factor  sinBi  in  x'  =  a'f  sinBj  describes  the  fact  that  two 
meridians  separated  by  a'  are  getting  closer  together  near  the  zenith  8=0 
(see  figure  5). 


The  quality  of  the  telescope's  lens  system  is  described  by  its  point 
spread  function  M  or  by  its  transfer  function  M: 


M(x',  y') 


e27ri(x'v+y'p) 


d  . 


(10.5) 


Hence  a  point  on  the  sky  at  ai  +  a1,  Bi  +  8*  will  produce  as  radiation  distri¬ 
bution  in  the  image  plane  the  point  spread  function  centered  at  x'  =  a'f  sinBj , 
y'  =  g'f.  In  mathematical  symbols, 


<$  (a  -  aj  -  a',  8  ~  Bi  -  8')  ->M(x'  -  a'f  singj ,  y'  -  8'f)-  (10.6) 
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Figure  5 
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If  the  point  source  gets  stronger  by  a  factor  of  say  SmjCoti  +  a',  +  3’)  then 

also  the  radiation  in  the  image  plane  gets  stronger  by  the  same  factor: 


All 

and 


smj(al  +  a'»  81  +  8')  6 (a  -  ax  -  a’,  B  -  3i  —  3') 

smj(ai  +  a'»  8i  +  8')  M(x '  -  a'f  sinBi,  y'-  B'f).  (10.7) 

these  point  sources  together  form  the  signal  on  the  sky,  whereby  the  spectral 
density  filters  Fm  and  Tj  are  already  taken  into  account: 


0!  +  3’)  6 (a 


“1 


a’,  3  -  3i  -  3')  da'dB' 
=  Smj(a,  3). 


(10.8) 


Because  the  telescope  system  is  linear  in  radiative  power  the  effect  of  Sm  j  (a ,  31) 
in  the  image  plane  is  computed  by  the  convolution  of  Smj  and  M: 


mj<al  +  a 


f 

* 


3i  +  3')  M(x ' 


-  a’f  sin3i ,  y'  -  3’f)  da'dB'. 

(10.9) 


Only  what  falls  onto  the  receiver  can  contribute  to  the  photoelectronic 
signal.  This  can  be  described  in  mathematical  terms  by  multiplying  the  above 
expression  by  the  receiver  weighting  function  R(x',  y')  and  then  integrating  it 
over  x',  y'.  The  R  is  "1"  where  the  receiver  is,  and  "0"  elsewhere.  For  the 
time  being  let  us  consider  only  one  receiver: 


+  a',...)  M(x'  -  a'f  sin3i,...) 

R(x',  y')  da'  dB'  dx'dy'. 


(10.10) 


The  terms  M  and  R  together  describe  the  angular  spread  function  P  of  the 
scanning  telescope: 


J J Smj(al  +  a'>  8l  +  8')  P (-a '  sinBi,  -3')  da'dB'; 

P(-a'  sinBi,  -3')  ~  I  I  M(x’  "  a’f  sinBi,  y'  -  3'f)  R(x',  y')  dx'dy' 

.  jf.  . 


.  .  27Ti(va'f  sin8i+y3'f)  ,  , 

-u)  R(v,  y)  e  dvdy. 


(10.11) 


Due  to  the  torsion  oscillation  the  azimuth  angle  ai  of  the  telescope 
axis  varies  according  to 

cq  =  a^  sin{27iJ2(t  -  to)}.  (10.12)  - 

During  time  intervals  long  compared  to  the  chopper  period  l/vc  we  can  approximate 
ai  ph  w(t  -  tg)-  (10.13) 

But  a)  will  be  different  in  different  time  intervals.  It  can  be  computed 
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as  the  derivative  of  ai(t)5 


to  =  a>(t)  =  do^/dt  =  2ttS7  aM  cos{27rfi(t  -  tQ) 


(10.14) 


Now  we  have  followed  the  signal  from  the  sky  into  the  receiver  where  it 


is 


2  Ai  cos (2ttvc t ) 


// 


^mj t  tg }  +  4  ,  8j  +  3*) 

P(-a'  sing!,  -6')  da’d3'  =  Umj(t).  (10.15) 


The  next  step  is  the  addition  of  receiver  noise  n(t),  which  is  supposed  to  be 
independent  of  the  signal,  and  specified  only  by  its  power  spectrum  |n(vt)|2. 
Signal  and  noise  together  are  subjected  to  an  electrical  filter  with  transfer 
function  G(vt)  before  being  recorded  on  tape  as  Up%(t) : 


uR(t) 


tumj  (O  +  n(t)  ]*  G(t)  = 

r  _ 

I  [Omi(vt)  +  G(vt)  e  TriVtt  (jVt> 


(10.16) 


§11  Influence  of  the  Neutral  Density  Filters  Tj 

The  neutral  density  filters  may  have  transmittance  factors  of  Tq  =  1, 

Tj  =  10-1  ,  T2  =  10~2 ,  for  example.  The  recorded  signal  will  suddenly  change  by 
a  factor  of  10  whenever  the  servo  demands  a  filter  change.  These  jumps  are 
disquieting  at  first  sight  for  two  reasons:  the  jumps  need  a  finite  amount  of 
time  during  which  we  loose  the  genuine  signal;  we  don’t  know  how  to  define  the 
average  signal  power  spectrum,  which  is  an  important  quantity  for  the  digital 
signal  restoration  process  (see  §5) . 

To  answer  the  latter  question  first  the  reason  for  switching  the  neutral 
density  filters  Tj  at  proper  occasions  is  of  course  to  maintain  a  suitable 
power  level  arriving  at  the  receiver.  In  other  words  the  signal  power  spectrum 
at  the  receiver  is  essentially  always  the  same,  no  matter  which  neutral  density 
filter  Tj  happens  to  be  in  operation.  The  short  times  it  takes  to  exchange  the 
Tj  filters  must  be  excluded  when  computing  the  average  signal  power  spectrum 
because  this  transition  has  nothing  to  do  with  the  genuine  signal. 

During  the  short  filter  exchange  times  the  signal  is  lost.  But  this  is 
not  severe  because  there  are  actually  four  receivers  in  the  image  plane  of  the 
telescope.  A  given  point  on  the  sky  is  observed  at  slightly  different  times  by 
the  four  receivers.  The  Tj  -  filter  may  or  may  not  be  switched  simultaneously 
for  the  four  receivers  depending  on  whether  the  pivot  point  of  the  filter  wheel 
is  in-line  with  the  four  receivers  or  far  off  the  line  of  receivers  (see  figure 
6) .  Even  more  generality  would  be  achieved  if  the  edges  of  the  filter  sections 
were  not  straight  lines  coming  from  the  pivot  point. 

In  any  event  when  plotting  in  an  a,t-diagram  the  dead  times  and  the 
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eadtimes  sequentially,  the 
a test  receivers  at  first 
his  is  the  best  case. 


resulting  dead  angles  for  the  four  receivers  (figure  6)  it  becomes  apparent  that 
it  is  preferable  to  let  the  dead  times  for  the  four  receivers  occur  sequentially. 
The  latest  receiver  should  be  switched  -first.  "Latest"  means  that  this  receiver 
scans  a  certain  sky  point  always  later  than  the  other  three  receivers.  However, 
the  property  of  "latest"  is  converted  to  "earliest"  after  the  torsion  oscillation 
has  reached  the  reversal  angle  a  =  ±  a^.  Furthermore  the  filter  wheel  has  to  be 
turned  sometimes  from  Tj  to  Tj+i,  but  also  sometimes  in  reverse  from  Tj+^  to  Tj . 
Hence  as  a  compromise  it  is  sensible  to  apply  simultaneous  switching.  The 
switching  dead  time  should  be  shorter  than  the  time  it  takes  to  scan  from  a  to 
a  +  Aa,  where  Aa  is  the  angular  separation  of  two  adjacent  receivers. 


§12  Dilute  6-Scanning 


The  available  observation  time  is  net  very  long,  allowing  onl’7  a  small 
number  of  a-revolutions ,  maybe  10  or  20.  Hence  the  elevation  steps  AB  must  be 
relatively  coarse,  say  5°.  This  is  much  more  than  the  diameter  of  the  angular 
spread  function  P  of  the  telescope  (lens  and  receiver  considered) .  This  has 
two  consequences:  we  don't  know  anything  abou.  object  points  midway  between 

two  6-scans ,  say  at  Si  and  at  62  =  61  +  AS.  Furthermore  our  measurements  are 
somewhat  ambiguous.  This  can  be  understood  with  the  help  of  figure  7,  where 
parts  of  /  P (a ' ,  6  -  Si)  da'  and  I  P(a',  6  -  62)  da'  are  plotted.  These  plots 
are  in  effect  weighting  functions  indicating  if  and  how  strong  a  sky  point 
a,  S  will  occur  in  the  recorded  signal.  The  ambiguity  becomes  obvious  when 
considering  two  object  points  at  019 »  61  and  at  ap,  61  +  e,  the  first  point  being 
weaker  than  the  second  one.  Nevertheless  the  recorded  signals  could  be  indis¬ 
tinguishable.  Even  a  line  object  of  the  form  6 (a  -  a0)  would  produce  the  same 
type  of  a  recorded  signal.  L.  Biberman  in  his  book  on  "Reticles  in  EJectro- 
Optical  Devices"  (Pergamon  1966)  also  discusses  this  ambiguity  in  length.  He 
concurs  that  this  ambiguity  cannot  be  removed,  except  by  decreasing  the  scanning 
step  AS  or  by  providing  some  additional  information  about  the  object  (see  §14). 


a 


Figure  7 
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iDfnw., .  ''""arrWTTflWTT  "WW 


Since  these  openings  don’t  exist  in  general  we  have  tc  admit  the 
ambiguity  and  try  to  make  the  best  of  it.  This  implies  that  we  might  replace 
the  unknown  signal  Sm . (04  +  a',  Si  +  S’)  by  an  equivalent  signal  SE(ai  +  a’, 
Si  +  S')  which  appears  to  be  the  same  to  cur  scanning  telescope: 


sinSi , 


-S’)  da'dS' 


S’)  da'dS*. 
(12.1) 


There  exists  an  infinite  variety  of  Sg  models,  from  which  we  may  choose  the  one 
which  is  most  convenient  in  terms  of  mathematical  manipulations  and  experimental 
test  procedures.  Two  of  the  most  obvious  choices  are 


Sm j  (a >  Si  +  3')  =  SD(a)  S(S’) 

Smj(a’  61  +  =  SE^a)  rect  (3 f /53) . 


(12.2) 


The  second  model  turns  out  to  be  very  convenient.  Herein  the  rect-function  is 
"one"  within  —6 S / 2  1  S’  <  6S/2  and  "zero"  elsewhere.  For  the  second  model  we 
get 

B  r  r  rr 

Smj  P  da'dS'  =  j  SE(ai  +  a')  rect(S'/63)  P(-a’  sinSi,  -S')  da’dS' 

SE(ai  +  a')  ,-ect(S'/6S)  P(-a'  sinSi,  -3')  dS'J  da’.  (12.3) 

In  the  square  bracket  we  can  disregard  the  rect-function  because  P  itself  is 
anvway  zero  outside  of  | 3 ' I  £  66/ 2.  The  square— bracket  is  therefore  the 
angular  line-spread  function  L: 


P (-a '  sinSi,  -S')  dS’  =  L(-o'  sinSi) 


(12.4) 


The  line— spread  function  is  easier  to  measure  than  the  point— spread  function 
because  test  objectsmay  be  infinitely  extended  in  the  S-direction.  L  can  be 
deduced  from  the  lens  transfer  function  M  and  from  the  spatial  spectrum  R  of 
the  receiver, 

ill-,  .  v  2ri(va'f  sin6i+y3’f) 

L(-a  sin3i)  =  1  /  jM(-v,  -y)  R(v,  y)  e 

JJJ  dvdydS  ’ 

1  (  S/  m  2riva’f  sinSi, 

=  —  J  M(-v ,  0)  R(v,  0)  e  idv; 


L(v)  =  J  L(a)  e"27T1Va  da  =  (f)"2  M(v/f,  0)  R(-v/f,  0) . 


(12.5) 


With  the  equivalent  object  S£  the  signal  (t)  in  the  receiver  (see  §2)  is: 


U  .  (t)  =  2  A,  cos  (2irv 
mj  1  1 


t)  ^ S£(o){t  -  t 


Q}  +  a’)  L(-a’  sinS-^)  da’.  (12.6) 
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Roughly  speaking  Umj  resulted  from  Sg  after  it  had  been  filtered  by  L  and 
modulated  with  chopper  frequency  vc.  The  next  steps  (12.7)  are  addition  of 
noise  n(t),  and  filtering  in  the  electronic  system  by  G(v^.) ,  and  recording  on 
tape  as  Uft(t).  These  steps  are  best  described  in  the  temporal  frequency  domain 


\(vc)  -  {6m.j(vt)  +  n (v t ) }  G(v£); 

6  .(v,)  -  /  U  ,(t)  e'2,rlvtt  dt  = 
mj  t  I  mj 


-A  II  +2iriv  t  c  t  -2TTiv  t  ,  ,  , 

Aj  I  I  e  c  L  e  t  da  dt  + 


w 


.)  £(Vt  ~  M 

/  \w  sinBj/ 


'v  -  v  \  ,  /v.  - 

t  c 


+2TTiV  to 


V 


+  S, 


+  V  \  /V  +  V  \  .  ~1 


(12.7) 


§13  Signal  Restoration  by  Digital  Filtering 


After  having  recorded  Ur  =  [Umj  +  n]*G  we  want  to  recover  from  it  as 
well  as  possible  the  original  signal  Sg(a).  The  term  "as  well  as  possible"  is 
meant  in  the  root-mean  square  sense.  The  recovery  procedure  will  consist  of  a 
linear  filtering  process  which  may  be  performed  digitally.  Before  actually 
determining  which  filter  function  D  is  optimal  we  will  recall  briefly  the 
classical  Wiener-filter  problem  and  generalize  it  somewhat. 

Suppose  an  object  signal  U(x)  has  been  corrupted  by  additive  noise  n(x). 
What  is  the  best  convolution  function  H(x)  for  object  restoration? 

f 

I  H(x  -  x’)  {U(x’)  +  n(x')}  dx’  «  U(x).  (13.1) 

Wiener  found  an  H(x)  which  achieves  the  job  optimal  in  the  least-mean-square 
sense.  For  computing  H(x)  or  its  Fourier  transform  the  so-called  Wiener-filter 
H(v)  we  must  know  the  power  spectra  |U(v)|2  and  |ft(v)|2  of  object  and  noise, 
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2 


x’)  (U(x’)  +  n(x’)}  dx’ 


dx  ->  MINIMUM; 


I  U  ( v) j  2 

H(v)  -  — - ; 

|u(v)  |2  +  | n( v) |2 


2tt1vx  , 
e  dv . 


(13.2) 


The  minimum  requirement  above  can  be  translated  into  an  equivalent  minimum 
requirement  in  the  Fourier  domain, 


U(v)  -  H(v)  (U(v)  +  n(v) |2 


dv 


MINIMUM. 


(13.3) 


For  this  "Fourier  translation"  we  made  use  of  the  Plancherel  theorem  and  of 
the  convolution  theorem 


2 


dv ; 


H(v)  V( .') 


27:ivx  , 
e  dv. 


(13.4) 


The  Fourier  version  of  the  minimum  requirement  is  very  convenient  for  a 
generalization  which  we  study  now. 

Let  us  assume  that  u  is  filtered  by  L  before  noise  is  added: 

u(v)  -*■  u(v)  L(v)  -*■  u(v)  L(v)  +  n(v).  (13.5) 


Now  the  modulus  of  the  minimization  integral  is 

|  u  ( v)  -  Hj(v){u(v)  L.  ( \»)  +  n  ( v ) }  |  *  |u(v)  -  Hj(v)  L(v){u(v)  +  n(v)/L(v)}|. 

(\3.6) 

Hence  the  Wiener-type  derivation,  which  led  to  H(v),  can  be  performed  as  before, 
except  that  h  -*■  h/L  and  H  -*■  fijL. 

(13.7) 

Now  let  us  be  one  more  step  'rai 1  additional  filter  process  G(v)  may  be 
part  of  the  system  after  the  ».  n: 

0  -*■  QL  -+■  QL  +  fi  -*•  OLG  +  f\G .  (13.8) 


1 0  |  2 


HjL  = 


a |2  +  | n/ 


Hi(v) 


|L(v)  u(v)|2 
L(v) { | Lu | 2  +  | n| ? ) 


We  look  now  for  a  filter  H2 ,  which  minimizes 

0  ( v )  -  H  -  (v)  G(v){n(v)  L(v)  +  h(v)}|2  dv  =  /  |0  -  H2G  L{fl  -  fi/L.)  |  2  dv. 

J  (13.9) 

This  has  again  the  "Wiener  shape",  but  with  ft/L  instead  of  n,  and  with  fl2GL 


/ 


* 
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u  2 


uL 


u 


+  n/L 


uL  2  +  n  2 


|u(v)  L(v)  |  2 


instead  of  H, 

H2G  L  = 

Ho  ^  v')  =  - - -  *  - — - 

G(v)  L(v)  |u(v)  L(v) | 2  +  | n(v) | 2  (13.10) 

(inverse  filter)  (Wiener  weighting  function) 

The  last  Wiener-generalization  has  already  almost  the  form  we  need  as  can 
be  seen  at  the  end  of  §12.  The  corresponding  terms  are:  SE  ~  U ;  L  ~  L;  G  ~  G . 
However  the  signal  term  SE  occurs  with  a  frequency  shift  of  ±  vc  due  to  the 
chopper  modulation.  Hence  we  first  must  demodulate  our  tape  recorded  signal. 

We  do  this  by  first  multiplying  it  with  2  cos{2tt  vc(t  -  t0)  and  then  rejecting 
by  means  of  a  low-pass  filter  all  terms  which  still  contain  vc  in  some  form. 


For  a  general  function  f(t)  -  /  f(v)  e 
it  , 2  cos{2tt  vc(t  -  t0)  means 


2irivt 


dv  the  process  of  multiplying 


~  ~  .  -2nitoVc  ,  2/  ,  \  +27titovc 

f(v)  -*•  f(v  -  vc)  e  u  c  +  f(v  +  Vc)  e 


(13.11) 


For  the  specific  function  of  interest  the  recorded  signal  {Um_.  +  n}*G  =  UR  means 
UR(vt)  =  {tJmj  (Vt-)  +  n(vt)  G(vt)  + 

Umj  (vt  “  vc)  G(vt  “  vc)  e  27Tit°Vc  +  Umj(vt  +  vc)  G(vt  +  vc)  e  °  °  + 

ft(vt  -  VC)  G(vt  -  vc)  e"27Tit°Vc  +  n(vt  +  vc)  G(vt  +  vc)  e  2^ 

We  insert  the  expression  for  Umj(vt)  from  §12 and  neglect  all  terms  with  ±  2  vc 
which  will  be  rejected  by  the  lowpass  filter 


Ai  e 


-27iivttr 


umj(vt)  =  u  sinBi 

_  /^t  “  vc\  „/vt  -  vc\ 
[•••]  "  sEy  ”  j  singly 


+2Trivct0  +  (same,  but  vc  -»•  -vc)  . 


(13.13) 

The  spectrum  of  the  recorded  signal  after  digital  demodulation  contains  as 
contribution  from  the  object 


or  abbreviated 
also  be  renamed 


(VJ\  -/  \ 

-\u  /  L\co  sinBi / 

*  se Ln  ^(vt)- 

d  for  brevity 


Al 


u>  sinBi 


e-27,ivtto  [6(Vt  _  V(;)  e-2”i'>cto  + 


G(vt  +  vc)  .  (13.14) 

The  noise  spectrum  of  the  recorded  signal  will 

-2nitnv, 

-  i-  n^,vt  t  vc; 

(13.15) 


brevity 

nE(vt)  =  n(vt  -  vc)  G(vt  -  vc)  e“2TTlt:oVc  +  n(vt  +  vc) 

-  .  +2iTitL  vc 

G (v^  "F  Vq )  e  . 
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It  is  important  that  the  effective  noise  power  spectrum  | fiE ( vt)  | 2  can  be  com¬ 
puted  from  the  original  noise  power  spectrum  | n ( | ^  and  from  other  known 
properties  such  as  the  electrical  filter  function  G(vt): 

lMvt>l2  "  -  vc)|2  |G(vt  -  vc)|2  +  |n(vt  +  vc)  | 2  |G(vt  +  vc)|2 

+  n(vt  -  vc)  n*(vt  +  vc)  G*(vt  +  vc)  e  ^7T^t0vc  4.  conj .  compl.  (13.16) 

Obviously  the  first  two  terms  are  known.  The  last  two  terms  can  be  considered 
to  be  zero,  at  least  whenever  a  vt-integration  takes  place.  The  phases  of 
h(vt  -  vc)  and  of  n*(vt  +  vc)  can  be  assumed  to  fluctuate  so  wildly  that  any 
integral  will  cancel,  or  at  least  fluctuate  itself  so  rapidly  that  the  demodu¬ 
lation  lowpass  filter  for  <  vc/2  will  block  any  contribution. 


At  this  point  we  found  the  cemporal  frequency  spectrum  of  the  recorded 
signal  after  demodulation  was 


UD(vt)  -  h{jr)  E(vt)  +  nE(vt).  (13.17) 

Herein  E  and  (figl2  are  known  and  Sg  or  its  Fourier  transform  Sg(a)  is  wanted  as 
accurately  as  possible-  The  modified  Wiener  filter  suitable  for  this  job  is 


SE(vt/w) | 2 

D(vt)  -  ■ST7~r  - -  . 

C  |sE(vt/m) | 2  +  |nE(vt)|2 


(13.18) 


§14  A  Method  for  Closing  the  Gaps  of  the  Dilute  6-Scanning 

In  §12  we  had  discussed  the  dilemma  which  was  caused  by  the  fact  that 
the  steps  A6  in  elevation  had  to  be  considerably  larger  than  the  resolution 
66  in  elevation.  The  factor  A6/66  might  be  about  five.  The  reason  for  the 
dilute  6-scanning  is  lack  of  observation  time.  Consequently  the  "image"  will 
have  wide  gaps  in  the  6-domain.  In  general  these  gaps  are  unavoidable  under 
the  given  constraints  due  to  the  scanning  mode  of  the  instruments.  However, 
there  might  be  special  circumstances  which  allow  us  to  close  the  gaps. 

Let  us  assume  there  are  only  point-objeots  on  the  sky.  We  propose  to 
replace  the  four  independent  receivers  by  a  single  receiver  which  might  have 
for  example  the  shape  of  a  V  or  an  N  (figure  8).  The  height  of  the  V  (or  the 
N)  Is  A6.  Hence  every  point  on  the  sky  will  have  at  some  time  the  opportunity 
to  send  some  light  into  the  extended  receiver.  We  will  know  of  course  the  A6- 
interval  for  every  object  point  since  we  keep  track  of  the  A6-steps  of  the 
scanning  procedure.  The  fine  resolution  within  a  given  6-interval  can  be 
deduced  from  the  fine  structure  of  the  recorded  signal.  For  example  a  point 
on  the  sky  which  is  scanned  by  the  lower  part  of  the  V  receiver  will  respond 
with  two  peaks  close  together,  while  another  point  a  little  higher  up  will 
produce  two  peaks  somewhat  farther  apart. 

The  mathematical  procedure  for  getting  fine  6-resolution  out  of  the  raw 
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signal  Sm(°0 ,  which  has  been  collected  by  scanning  with  the  V-receiver  over 

the  strip  (3m>  3m  +  A3)  is  as  follows:  Suppose  we  want  to  compute 

S(ct,  3m  +  i<53) ;  i  =  0,  1,...,  I  —  1 ;  163  =  A3»  The  shape  of  the  V  may  be 

such  that  the  two  bars  are  apart  from  each  other  by  (i  +  2) 6a  at  an  elevation 

i3  above  3m.  If  there  is  a  point  object  on  the  sky  at  (a,  3m  +  i<$3)  then  it 

will  create  two  pulses  in  the  raw  signal  Sm(a)  at  a  +  (i  +  l)6a  and  at 

a  -  (i  +  l)6a.  We  test  this  by  multiplying  the  signals  Sm{a  +  (i  +  l)6a}  and 

Smiot  (i  +  1) 6a} .  For  adjusting  the  magnitude  of  this  object  point  we  have 

to  draw  the  root 

SSm{a  +  (i  +  l)6a}  Sm{a  -  (i  +  l)6aj  =  S(a,  Sm  +  i63); 


i  =  0,  1. . . ,  I  -  1;  163  =  A3. 


(14.1) 


In  order  to  achieve  3~resolution  in  63  steps  we  have  to  repeat  this  operation 
for  all  possible  double  pulse  separations.  Now  S(a,  3m  +  i63)  is  ready  for  the 
filtering  process  described  in  §13. 


Possibly  this  filtering  operation  might  have  to  be  generalized  into  two 
dimensions.  But  such  a  study  is  not  warranted  at  present  since  it  is  unknown 
if  all  or  most  of  the  objects  are  points,  and  if  there  exists  any  chance  to 
implement  a  V-receiver.  However,  if  these  two  "ifs"  are  confirmed  an  additional 
study  is  desirable  for  finding  what  happens  if  there  are  occasionally  objects 
which  are  not  points . 
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V' 

It  is  planned  to  send  an  infrared  scanner  in  a  balloon  up  to 
high  altitudes.  Azimuth  scanning  is  accomplished  by  the  torsion  of 
the  gondola.  The  elevation  angle  of  the  scanning  telescope  is  incre¬ 
mented  stepwipe  after  each  full  revolution.  In  this  report  the 
scanning  mode  is  studied  and  improvements  thereof  are  suggested. 

For  a  follow-up  instrument  it  is  suggested  that  the  scanning  rate 
should  not  be  limited  by  the  slow  torsion.  In  a  next  approach  the 
elevation  would  be  scanned  fast  compared  to  the  torsion.  Ultimately 
a  "Multiplex  Scanner"  should  look  at  many  directions  simultaneously. 
The  principle  of  a  multiplex  scanner  is  outlined. p 


.V  I  ~  r  fHff.  a  r 


